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Abstract 

The Singlet-Doublet model of dark matter is a minimal extension of the Stan¬ 
dard Model with dark matter that is a mixture of a singlet and a non-chiral 
pair of electroweak doublet fermions. The stability of dark matter is ensured 
by the typical parity symmetry, and, similar to a ^Bino-Higgsino’ system, the 
extra matter content improves gauge coupling unification. We revisit the ex¬ 
perimental constraints on the Singlet-Doublet dark matter model, combining 
the most relevant bounds from direct (spin independent and spin dependent) 
and indirect searches. We show that such comprehensive analysis sets strong 
constraints on a large part of the 4-dimensional parameter space, closing the 
notorious Tlind-spots’ of spin independent direct searches. Our results empha¬ 
sise the complementarity of direct and indirect searches in probing dark matter 
models in diverse mass scale regimes. We also discuss the LHC bounds on such 
scenario, which play a relevant role in the low mass region of the dark matter 
candidate. 
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1 Introduction 

In defiance of a plethora of experimental searches, the mass and dynamics of Dark Matter 
(DM) remain largely unknown. Up to now, there is no indication that the DM particle 
interacts with the known matter other than gravitationally. However, if one assumes 
that the relic abundance of DM today is determined by the same physical processes that 
determined the relic abundance of baryonic matter, one would require that DM at some 
early point in the history of the universe was in thermal equilibrium with the Standard 
Model (SM) particles. 

The idea of a thermal relic DM had been further motivated by the following observation: 
A particle whose mass and dynamics are dictated by electroweak (EW) physics would have 
a thermal evolution that can reproduce the observed DM relic density today. In simple 
constructions of this type, the DM has spin independent interactions with the nucleons via 
a Z-boson exchange that are orders of magnitude larger than experimental limits. Another 
popular proposal that relates the mass and dynamics of DM with the electroweak sector 
is the neutralino of supersymmetric extensions of the SM. There, the mass of the DM 
candidate is related to the supersymmetry breaking scale, which in turn is related to the 
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EW scale by means of the naturalness criterion. In this case the DM is a Majorana particle 
so that the spin independent DM-nucleon interaction via Z-boson exchange is absent. 

The negative searches in DM and collider experiments for these explicit scenaria has 
led to the consideration of different approaches in describing the properties of DM. One 
direction that has been extensively studied is by using the tools of effective field theory, 
where the coupling of DM with SM particles is parametrised by higher dimensional opera¬ 
tors suppressed by the scale that determines the validity of the effective theory. While this 
approach has been very useful in obtaining model independent bounds on DM properties 
from direct and indirect searches, it has been shown mum that it can be misleading 
when considering the complementarity of these searches with DM searches in high-energy 
colliders. 

In view of this difficulty, simplified Models of DM have been lately motivated as an 
alternative to effective DM theories that allow for a unified study of the multiple DM ex¬ 
perimental fronts pirns]- Furthermore, the negative LHC searches for complete ultraviolet 
extensions of the SM with a natural DM candidate (e.g. supersymmetry), point towards 
model building beyond the schemes that specific UV completions impose. In this perspec¬ 
tive, a well motivated simplified Model of DM can capture the dynamics of different UV 
completions that are relevant for each specific experimental search. 

In this paper we discuss the theory and phenomenology of a Simplified Model of DM 
that consists of a singlet and an electroweak doublet fermion mum, besides the SM field 
content. The dark matter candidate of the ‘Singlet-Doublet’ Model (SDM) is identified as 
the lightest neutral fermion of the new sector whose stability is ensured by a Z 2 symmetry. 

Within the vast landscape of DM models, the SDM enjoys certain motivating features: 
i) Given the tight experimental constraints on chiral dark matter, it represents a minimal 
model of electroweak dark matter that is stable by a symmetry of the theory]^ ii) The 
sensitivity of current direct search experiments is high enough that DM-nucleon interactions 
via a Higgs boson exchange are now being probed. The SDM is a simplified model that can 
be used to parametrise these bounds since the DM candidate has tree-level interactions with 
the Higgs boson, iii) The new field content of this model improves the unification of gauge 
couplings with respect to the Standard Model (however the unification is reached at a scale 
which is too low to avoid proton decay. A way to address this issue has been proposed in 0 )- 
iv) The SDM captures the relevant dynamics of well-studied supersymmetric candidates 
such as the Bino-Higgsino system in MSSM and the Singlino-Higgsino system in NMSSM- 
type scenarios. 

^For a similar construction with an electroweak doublet and a triplet, see [5]. 
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Many aspects of the theory and phenomenology of SDM have been investigated in the 
past [3 na im [El uni na nn], including the recent analysis of [H] where the constraints 
on the parameter space have been updated focusing on spin independent direct detection 
experiments. 

The goal of the present work is to explore the complementarity of dark matter and 
collider searches for the SDM, especially in view of the vast progress that has been recently 
achieved in the indirect search observational front. In particular, we see how a combination 
of spin dependent and indirect searches probes regions of the parameter space where spin 
independent searches lose sensitivity, the so called ‘blind spots’. Combining these bounds, 
one can reach sensitivities that allow for exclusion of a very large part of the parameter 
space, extending up to DM mass values of the order TeV. 

We also focus on the special features that the model exhibits when the DM candidate 
is relatively light, in particular on the extra constraints coming from EW precision observ¬ 
ables as well as Z-boson and Higgs boson decay width. In particular the LHC constraint 
on Higgs invisible decays has a relevant impact on previously unbounded regions of the 
parameter space. We will observe that the low mass scenarios with relatively large Higgs- 
DM interaction are almost completely ruled out, except for small corners of the parameter 
space that will be accessible at the next LHC run. 

The structure of the paper is the following: In section we introduce the model, we 
discuss the relevant DM couplings to SM particles and we explain the main DM annihilation 
channels. In section we display the limits set on the parameter space of the model by 
indirect and direct DM searches, showing their complementarity. In section]^ we focus on 
the low mass region and show the interplay of DM searches and collider constraints and 
finally, in section we summarise our results. 

2 The Model 

We now introduce the Singlet-Doublet Dark Matter model (SDM) that we will study in 
the rest of the paper. The field content of the dark sector consists of new Weyl fermions 
with the following transformation properties under U(l) and SU(2): 

A: (1.0). : (2.-1). V-2 = : (2. 1) (1) 

The electroweak doublets have opposite hypercharge in order to cancel gauge anomalies as 
well as to avoid the strong bounds on chiral dark matter [IS] . Furthermore, the new fields 
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are assumed to be odd under a Z 2 symmetry under which all SM fields are even. With 
this field content, apart from the obvious couplings to gauge bosons from the covariant 
derivative, interaction with the SM via renormalisable operators is achieved exclusively by 
coupling to the Higgs. The dynamics of the new fields are then governed by the following 
Lagrangian: 

_ _ _ JYl 

J^SDM = 2 ( 

-yiiix-HX-y 2 'ip 2 ^H + h.c., (2) 

and determined by the four parameters 

ms, mn, ?/i, 2/2, (3) 

denoting respectively the mass terms of the singlet A, the doublets and the Yukawa 
couplings among these fermions and the SM Higgs doublet H which can also be expressed 
as 

2/1 = 2 /cos 6», 2/2 = 2 /sin 6*. (4) 

In our analysis, we will neglect the overall CP violating phase of these parameters and take 
them real, while, without loss of generality we also take mg, m^, y > 0. For constraints 
on SDM from limits on the electron EDM dg in the case of a non-vanishing CP violating 
phase, see [0]j^ 

Upon electroweak symmetry breaking, the neutral components of the doublets and the 
singlet mix. The physical spectrum of the model then consists of three Majorana fermions 
with increasing mass mi^ 2,3 and a Dirac charged fermion whose mass we denote as m^+. 
At tree level, one has m^+ = —m^. The mass matrix of the neutral states and the rotation 
to the mass eigenstates are given by 
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(5) 


The lightest eigenstate xi is our dark matter candidate]^ Its composition in terms of the 

^Considering the recent updated limit \de\ < 8.7 x 10“^^ ecm [17], we do not expect the improved limit 
on the single physical phase (that can be chosen to be that of mo) to be below 0(0.1) even for ms, mo of 
0(100) GeV. 

^We use X and xi (^s well as mi and m^) interchangeably. 
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fields in the gauge basis is 


X, = UnA + \Un? + \Uu? + \Ui3? = I- (6) 


The dark sector interacts with the SM via couplings to the W, the Z and the Higgs 
bosons. Written in the mass eigenstate basis and standard 4-component notation]^ they 
are given by 


C D hX Zi^X {cZxmXnPL CzxmXn^^)^'^ 

-T(K.3ir;V.7'‘n'i'+ - c/iw;V„ 7 '‘P„<i<+ + hx .), 


(7) 


where 


- U^2K2) , ChxmXn = ^{yiU*m2Unl + y2U*r^3Kl) • (§) 


The couplings of the DM candidate then take the simple form 


- Ut^W-X^^>^PR^+ + h.c .), 

where Chxx Czxx can be written as 

^ {‘iymm^ P {yl P yl)mi)v 

n^-D + ( 2/1 + 1 / 2 )T + 2?n-s ?n-i - 3mf ’ 

^_ mzvjyj - yl){ml - ml) _ 

2{m\ - m2 )2 + t)2 mi rrio + (yl + 2 / 2 )(”^i + f^D) ' 


(9) 


( 10 ) 

( 11 ) 


We note that mi here can have either sign. From these expressions we observe that there 
are corners of the parameter space where the above couplings vanish. Such “blind spots” 
of the Higgs PE! and Z boson couplings to dark matter occur for 


sin(26*)mD -|- mi = 0 

I2/1I = I2/2I or |mi| = mD 


Regarding the coupling to the Z boson, the blind spot at 


^We have = {xna x")^ and 4'+ = (V’l 


c/ixx ~ 9; (12) 

czxx = 0 • ( 13 ) 

1 2 / 1 1 = 1 2 / 2 1 is easily understood by 
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looking at of oq. (11): when the two Ynkawa couplings are equal, |17i2p = so that 

the coupling vanishes. This is analogous to the vanishing of the neutralino coupling to the Z 
at tan/3 = 1 in supersymmetric scenarios. The second possibility, |mi| = nio, corresponds 
to the limit where DM is almost purely made of the neutral components of the doublets, 
thus again to a case where the two mixing angles are almost equal |17i2p — — 1/2. 

In other words, in this second case c^xx vanishes as DM becomes a Dirac fermion and the 
vector interaction of standard Dirac DM is restored. 

Finally, we remind that the Lagrangian of eq. 0 represents a generalisation of the 
Bino-Higgsino sector of the MSSM in the decoupling limit of only one light SM-like Higgs 
boson. In the MSSM, supersymmetry dictates that the Yukawa couplings yi and y 2 are 
related to the U(l) gauge coupling g’. Instead, we take them here as free parameters. The 
MSSM limit of the model can be obtained as follows: 


y 


g' 


rriQ Ml, rUn —?• —/r, cos 9 ^ — cos/3, sin 9 sin/3, 


where Mi is the SUSY-breaking Bino mass term, /r is the Higgsino mass term and tan/3 
the ratio of the vevs of the two Higgs doublets. 


2.1 Dark Matter annihilation channels 

Before we move on with the discussion of the constraints to SDM, we give here a brief 
overview of the main annihilation modes of the DM candidate. The annihilation modes 
participate both in the determination of the thermal relic density in the early universe and 
in the expected indirect detection rate in the late universe. 

In the following, for each annihilation channel, we focus in particular on the correspond¬ 
ing contribution for indirect detection. We note that the annihilation of Xi to fermions in 
the late universe can only be mediated by a Z boson, since the s-wave amplitude of the 
Higgs-mediated contribution vanishes. Due to the helicity structure, the s-wave contribu¬ 
tion to the Z-mediated annihilation cross section scales with the fermion mass. Therefore, 
in summary, the main annihilation channels in the SDM are to bosons and for high enough 
DM mass, to tt. 

Dark matter annihilation into W~^W~ occurs at tree level via a t-channel exchange of the 
charged fermion and via an s-channel exchange of a Z or a Higgs boson. However, the 
former dominates the rate relevant for indirect detection since the Z/h exchange has no 
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s-wave contribution. Annihilation into ZZ occurs at tree level via a t-channel exchange of 
a neutral fermion Xi via an s-channel exchange of a Higgs boson. As before, only the 
t-channel is relevant for indirect detection since the s-channel does not contribute in the 
vanishing velocity limit. Sommerfeld enhancement brings corrections to the tree level rate 
in the limit ^ raw From the results of [18], we expect that the enhancement appears 
for Xi that is almost purely doublet and is marginal in the parameter space we consider. 
We will therefore for simplicity neglect this effect. 

XX ^ ti 

Dark matter annihilation to fermions occurs in an s-channel, through a Z or a Higgs boson. 
Because of the helicity structure of the process, the Higgs channel is irrelevant for indirect 
searches while the s-wave amplitude for the Z channel scales with the fermion mass so 
that, for heavy enough dark matter, the annihilation to two top-quarks is the dominant 
one. 

XX ^ Zh 

This process proceeds via t-channel exchange of a neutral fermion Xi via s-channel 
through a Z boson. Both channels can potentially be relevant for indirect searches since 
none of them is velocity suppressed. Since the t-channel passes through a XiXjh vertex, 
we expect that the rate of this process grows with y (unless it falls on a blind spot). 

XX ^ 77, iZ 

Annihilation into two photons or a photon and a Z boson occurs at one loop via box 
diagrams of W bosons and and via an s-channel Z and a triangle loop of SM fermions, 
or W bosons. The calculation of the loop induced annihilation rate of DM to 77 or 
7 Z in the MSSM context with a neutralino DM candidate was done in [T^ [20l [ 21 ]. It 
was shown that for Higgsino DM heavier than around 300 GeV, the cross section for both 
channels reaches a plateau at ~ 10 “^® cm^s“^ while the non-perturbative, non-relativistic 
corrections introduce a resonant peak at ~ 6 TeV [22] 

Other annihilation channels include xx ^ hh and xx //, where / are SM fermions 
other than the top. Since the DM particles are Majorana fermions, the first one is again 
velocity suppressed and thus insignificant for indirect searches. As mentioned above, the 

®In the parameter space that we consider in our numerical analysis in Section]^ (with TOd, nis < 5 TeV), 
the non-perturbative enhancement can be at most of one order of magnitude in the annihilation cross 
section (see Fig. 5 of jH]). 
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annihilation to fermions pairs, in particular bb and r~^T , is mediated through a Z boson 
and helicity suppressed, therefore possibly relevant only for DM masses below mw 

3 Limits from Direct and Indirect Searches 

In this section, we use a collection of direct and indirect searches for dark matter to put 
constraints on the parameter space of SDM. The bounds from spin-independent direct 
detection have been recently investigated in [8] where it was shown that, even if they are 
strong for relatively large Yukawa couplings, they are ineffective around the ‘blind spot’ 
regions, where the coupling of Xi fo Higgs vanishes. We will see that, after taking into 
account a combination of spin-dependent and indirect searches, these blind spot regions 
can be probed and excluded. 

Throughout our analysis, the thermally averaged annihilation cross section of xi info 
the channels discussed in the previous section have been computed by modifying appropri¬ 
ately the expressions that can be found for example in [2H] • We cross-checked our analytical 
calculations employing micrOMEGAs |2lj, finding an excellent agreement. 

In the following analysis, we assume that the lightest state xi accounts for the whole 
dark matter abundance. Nevertheless, we show in the plots a line where its relic density 

as computed by micrOMEGAs agrees with the Planck observations, — 0.12 

[25]. Such a line divides the parameter space in two regions, where xi as standard thermal 
relic would be respectively under and overabundant. Thus, xi can be considered a good 
DM candidate (and the bounds we are going to show apply) in the former case if a non- 
thermal production mechanism is at work in the early universe, while in the latter case a 
non-standard thermal history of the universe, providing DM dilution, has to be assumed. 
If, on the contrary, one sticks on xi as a standard thermal relic, the areas of our plots 
leading to overabundance would be excluded by the results of Planck on while in 

those with underabundance Xi would be interpreted as a subdominant component of DM 
and the sensitivity of DM experiments to our model would decrease due to the reduced Xi 
density. We will comment about this possibility when discussing our results. 

3.1 Experimental input 

In our analysis, we use recent results by the Fermi-LAT (22], PICO [22], IceCube [2H] and 
LUX [29] collaborations, including prospects for the CTA experiment [30]. Before we move 
on with the results, we briefly report a few details about these searches, starting from 
searches for gamma rays. 
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The gamma ray flux obtained from annihilation of DM is given by 

1 lrFAv\ r^rnax _^ P P 

$,(Ar7) = — Y,Bi^dE,x / pl^{i,n)didn, (i4) 
»7r Je^.^ . a^j Ji.o.s.JAn 

J 

where (xav) is the velocity averaged total DM annihilation cross section, is the DM mass 
and BidNydE^ is the differential gamma ray yield for a particular final state multiplied 
by the corresponding branching fraction. The J-factor encodes the information about 
the astrophysical density distribution of DM at the source. It is defined as the integral 
along the light of sight (l.o.s.) of the square of the density of the DM distribution in 
the source, integrated over the solid angle of observation. Under the assumption of a 
particular DM density distribution profile, the flux upper limits can be used to set upper 
bounds on (cr^u). Dwarf spheroidal galaxies have the advantage that the integrated J- 
factor is relatively insensitive to the choice of a DM density profile isnisg. The bounds 
on the annihilation cross sections of % are obtained from the recent results of Fermi-LAT 
[2S] where the negative search of gamma rays coming from 15 dwarf spheroidal satellite 
galaxies of the Milky Way is used to set upper limits on the DM annihilation cross section 
into certain channels, namely e+e“, , r~^T~, uu, bb and WW. 

Searches for antiprotons have shown to set limits on %% WW/ZZ (which are in 
most of the parameter space the dominant annihilation channels) that are either smaller 
or comparable to the ones set from 7 -rays, depending on astrophysical details of their 
propagation through the galaxy [33l [Ml [35]. We do not consider these searches in the 
present study, however, it is certainly interesting to explore to what extend these searches 
can improve the bounds presented here. 

Concerning direct DM searches, we employ bounds on the spin-independent (SI) y- 
nucleon scattering cross section and the spin-dependent (SD) y-proton scattering 
cross section coming from LUX [20], PICO 1271 and IceCube [2H|- The latter one is 
an indirect, model-dependent limit on obtained by interpreting the bounds on the 
flux of neutrinos from DM annihilation in the sun. Two different upper limits are given in 
|28j . a stronger one by assuming DM annihilation to WW (resulting to neutrinos with a 
harder spectrum) and a weaker one by assuming annihilation to bb (resulting to neutrinos 
with a softer spectrum). In our case, for > mw, apart from certain regions of the 
parameter space where DM is mostly singlet, the dominant annihilation leads to hard 
neutrino spectra. One can recast the IceCube bounds for multiple annihilation channels 
with less or more hard neutrino spectra. Instead, in this work we take a conservative 
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Figure 1: Limits and CTA 

prospects on SDM for y = 0.01 
and any tan0. The red region 
is excluded by Fermi-LAT limit 
[2S] on {av)ww- The prospective 
exclusion by CTA (from [30]) is de¬ 
noted by the dashed grey line. The 
excluded regions are drawn under 
the assumption that matches 
the observed value (see text for 
details). Regions where y wonld 
be overabundant (nnderabundant) 
are denoted by “over” (“under”). 


approach and apply these bounds only when the annihilation to hard spectra is dominant 
[SS], i.e.: {(jv)ww + {(^v)zz + > 10 

3.2 Results 

The complementarity between different searches is illustrated in the aggregated plots of 
figures m In these plots we show the combination of bounds from indirect and direct 
searches as well as the prospects for the CTA experiment for a wide range of mass param¬ 
eters, while the Yukawa couplings yi ^2 span from small valnes (we show a representative 
case with y = 10“^ but the constraints are the same for smaller values) and ‘MSSM-like’ 
values [y = 0.2) up to large values substantially enhancing doublet-singlet mixing effects 
(we choose y = 1, a valne still compatible with pertnrbativity up to the unification scale 
|6|.). As it is clear from the mass matrix of the model, the relative hierarchy of the Yukawa 
couplings does not have a physical significance so it is sufficient to restrict to | tan0| ^ 1. 

Regarding constraints from gamma rays, as we discussed above, among all the annihi¬ 
lation channels the strongest bounds come from the upper limits on the annihilation cross 
sections yy WW, ZZ, it and Zh. In the following we discuss each channel separately. 

The regions excluded by Fermi-LAT for the WW channel are depicted in pink in figs. 

As discussed in the previous section, the annihilations to WW and ZZ are driven 
by electroweak interactions, hence they only depend on the doublet component of y, see 
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eq. @ jj Thus the resulting bound does not diminish in the limit of small Yukawa couplings, 
where it gives the only constraint to the model as shown in fig. [^for y = 0.01. In general 
the limits from %% —?• WW constrain a region of the parameter space where a thermal 
relic would be under abundant. Thus y would be a viable DM candidate only under the 
assumption of an additional non-thermal production. If on the contrary we assume that 
X is a thermal relic and therefore a subdominant DM component, one should take into 
account the reduced y density and rescale the bounds by a factor (Q^h^/0.12)^, resulting 
in practically no constraints in almost all the underabundant part of the parameter space. 

In the summary plots of figs. [T]j^we have also shown the prospects for yy WW from 
the CTA experiment, as estimated in PI- It is interesting to note that the improvement 
is only marginal and in some regions of the parameter space (e.g. for < 200 GeV) the 
sensitivity seems to be comparable or even worse than that of Fermi-LAT. This happens be¬ 
cause the improvement by CTA is substantial mainly in the high DM mass region where the 
annihilation cross section to WW is anyway too small (unless it sits on a non-perturbative 
resonant peak). 

The annihilation to tt is also important once this channel becomes kinematically ac¬ 
cessible. The latest Fermi-LAT results [26] that we use do not report on limits on dark 
matter annihilating to ti. However, since the photon flux from tt follows closely the one 
of bb |3Z|, for the purpose of demonstrating the bounds we find it sufficient to use the 
limits on the latter. The corresponding bounds are shown as purple regions in figs. |^and 

The annihilation occurs via s-channel Z boson mediation and thus the amplitude is 
proportional to the coupling czxx ~ ~ I^i 2 p- As we discussed below eq. (13), the 

channel is suppressed when y is either almost a pure doublet or pure singlet. Indeed, we 
see in figs. and that for small values of yi ^2 where the mixing between the singlet and 
doublet neutral states is small, the annihilation to tt, even when kinematically allowed, is 
highly suppressed so that it does not give any bound. For larger y, as long as | tand| is not 
too close to 1, yy ^ ft sets substantial bounds on the parameter space (see purple regions 
in fig.||), again if y is a non-thermally produced DM candidate, as previously discussed for 
the WW channel. These bounds, along with ones coming from the SD searches, can probe 
the blind-spot regions of SI searches that are controlled by the coupling to the Higgs Chxx- 

The annihilation channel yy —)■ Zh has a relatively different behaviour than the 
WW, ZZ and tt channels in the sense that its size is bounded from above by the Yukawas 
instead of the electroweak coupling. Therefore, it plays an important role for large values 


®The photon fluxes from yy —>■ ZZ are similar to those from yy WW therefore the corresponding 
limits are similar. Since in our model the annihilation cross section to WW is over almost all the parameter 
space larger than the one to ZZ, we show limits from yy —)■ WW. 
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Figure 2: Limits and CTA prospects on SDM for y = 0.2 and various values of tan0. Red: 
region excluded by {(tv)ww [2S]- Blue: exclusion from direct detection limits on [22! • 
Magenta: exclusion from IceCube limits on |2H]- CTA sensitivity prospect is shown 
as a dashed grey line. The excluded regions are drawn under the assumption that 
matches the observed value (see text for details). Regions where y would be overabundant 
(underabundant) are denoted by “over” (“under”). 


of y. For such y the bound on SI scattering cross section from LUX becomes very strong. 
Nevertheless, the constrain from yy —)■ Zh is seen to close the blind spot regions of the 
LUX bound. In particular, for y = 2 and tan6* = —2 (not shown in the plots) it excludes 
the region 20 < mg/GeV < 400, 200 < mo/GeV < 440, part of which is not constrained 
by LUX. 

The latest results of (22] include strong bounds on DM annihilating to light fermions: 
For yy —>■ ,bb, the latest limits reach as low as ~ 5 x 10“^^ cm^ for < 10 GeV, 
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Figure 3: The same as before for r/ = 1. Additional excluded regions that appear are 
the following. Green: exclusion from direct detection limits on m Purple: region 
excluded by {o'v)^ [26] . 


much below the thermal relic cross section. However, due to the helicity suppression of dark 
matter annihilation to light fermion in our model, the cross section is generally smaller: 
{av)rT ^ 3.5 X 10“^® cm^ s“^, < 5 x 10“^^ cm^ s“^. For this reason there are currently 

no bounds from these channels. 

The annihilation to two photons or a photon and a Z boson is maximised for % 
that is mostly doublet. Similar to the supersymmetric case of a pure Higgsino DM, 
the thermally averaged annihilation cross section reaches a constant maximal value of 
{av)-y-y^Z'y ~ 10“^® cm^ s“^ in the limit of large while the full non perturbative calcula¬ 
tion includes resonances that can increase this value by several orders of magnitude for a 
small mass splitting between % and the charged state '4)'^ [THl [ElEOlHUES]- However, in 
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Figure 4: Same as in the previous figure in the plane {y, ttiq — rrij^) for various choices of 
tan^ and m = (rUg + mY))/2. The blue fines indicate contours of the DM mass m^. 


the parameter space on which we focus these effects are seen to be small enough to evade 
constraints on the model. 

In addition to annihilation of x, the limits on spin-independent and spin-dependent 
scattering cross sections set strong constraints on the model, as seen in the plots of figs. 

1^ The SI interaction is mediated at tree level by a Higgs boson and starts to give important 
constraints already for moderate ‘MSSM-like’ values of y, as shown in fig. In particular, 
we notice that, in contrast to indirect searches, bounds from direct detection are relevant 
also for thermal x with ~ 0.12, i.e. in the region with mg ~ mo where singlet-doublet 
mixing is enhanced |8]. In this case LUX can set limits on up to ~ 1 TeV. Furthermore, 
as seen in fig. for large y the bound from LUX becomes very strong, excluding x as heavy 
as several TeV in the region of enhanced interaction, where Uu ~ U \2 ^ U 13 . We also note 
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that this interaction exhibits “blind spots”, i.e. regions of the parameter space where Chxx 
becomes identically zero, as we discussed in section In such regions, as well as when 
y is small {y < 0.1) the spin independent searches pose no bound on the model. On the 
other hand, the SD interaction, even if it is relatively less constrained by observations, it 
does not exhibit the same blind spots, as it is mediated at tree level by a Z boson. Hence, 
similarly to %% —?• tt, is enhanced at large y as the mixing becomes bigger and broader 
in the parameter space. The bound set by PICO, shown in green, is especially relevant 
in the low mass region. Furthermore, as we mentioned above, for DM heavier than mw 
and primarily annihilating to final states that decay to hard neutrino spectra, the indirect 
bounds on spin-dependent cross section from IceCube can be employed and, as seen in 
figi they lead to strong bounds depicted as magenta regions. As previously mentioned, 
both these bounds on as well as the indirect one on (crn)^^ can cover certain blind 
spots of LUX. This can be seen for example for ta,n9 = —2 in fig. 

The coverage of the low-y and blind-spot regions is also illustrated in fig. where we 
see how these bounds change for different scales of y. The plots are made in the plane 
{y, m-s — mn) for different values of tan0 and the average singlet and doublet mass term 
m = {ms -|- 7710 ) 72 . In the top two plots we focus on light DM, taking m = 200 GeV, 
corresponding to mx in the range 100-r200 GeV. We observe that as y increases, the 
LUX blind spot becomes broader and for 7 / > 2, SI searches set no bound on most of the 


parameter space, in particular when tand is negative and rUg < ruo. As we see, this part is 
covered by SD direct and indirect searches from WWjZZ^ tt and IceCube. The coverage 
of blind spots by these searches is also seen in the bottom two plots, where DM can be as 
heavy as ~ 600 GeV. 

Putting all the complementary exclusions together allows one to set bounds on % for a 
wide range of the parameter space. Here we focus on > 100 GeV as we will treat the 
special features and details of the low mass region in a separate section below, including 
also constraints from colliders. The bounds for the benchmark points of y and tand can 
be directly extracted from the contours iu the plots of figs. [T]j^ Scauuiug over differeut 
values of ism 6 for y < 0.01 has no measurable effect so that oue cau exclude % with 


mx A 280 GeV aud Uh < 0.5 for auy value of tau6*. For the ‘MSSM-like’ value of 7 / = 0.2, 
the limit depeuds ou the iuterplay betweeu the SI aud SD bouuds as we vary taud, which 
iu turu is determiued by the depeudeuce of Chxx czxx ou taud, as cau be iuferred from 
eqs. (10) aud 0. For positive tau 6 the SI bouuds are stroug as it cau be seeu iu the right 
plots of fig. For taud < —1.5 aud as we decrease the valueQat first the domiuaut bouud 


^At the special point tan 9 = —1 the mixing of the dark matter candidate practically vanishes so the 
situation is similar to that of the small Yukawa coupling. 
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Figure 5: In this illustrative 
18 scatter plot of avww ver- 

16 sus for y = 0.01, we 

14 clearly see how the annihila- 

tion cross section to WW is 
g enhanced as the mass split- 

\ ting becomes smaller. The 
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^ blue line depicts the Fermi- 

® LAT limit [2S]- 
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0 
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is the SD (from IceCube) but at around tan 6 * = —4 it gives its place to the SI bound. All 
in all, for y = 0.2 we obtain the exclusion < 220 GeV for Uh < 0.65. For larger Yukawa 
couplings, this interplay between the various searches is even more apparent. In general 
we obtain stronger bounds: For r/ = 1 we obtain that y is excluded for ^ 275 GeV 
unless it is almost a pure singlet < 0.8). For positive tan 6 * this bound becomes much 
stronger: For example, y is excluded for < 1 TeV unless it is a purely singlet or purely 
doublet state (0.015 < < 0.95). 

Exclusion on y — mass degeneracy 

The more xi is doublet-like, the more it will have a mass that is nearly degenerate with 
X 2 and . Therefore, the constraints on the mixing angles that we discussed above can 
be directly translated into exclusions of the y — mass splitting. 

This is another set of constraints where indirect searches play a major role. Indeed, 
indirect searches can be used to close in on the parameter space of scenaria with compressed 
spectra, where other experiments such as colliders and direct searches lose sensitivity. 
In particular, the SI and SD cross sections are smaller when xi is almost degenerate 
with the heavier partners because the mass degeneracy is generally an outcome of Xi 
being mostly doublet, so that |Viip ~ 0 and |? 7 i 2 p — |Vi 3 p ~ 0.5. In this case, both 
the ~ ViU 12^11 + I/ 2 V 13 V 11 and the ~ |Vi 3 p — |17i2p couplings are small and 
therefore direct searches cannot probe the highly degenerate regime. On the other hand, 
the annihilation to WW/ZZ is enhanced in this regime since these channels are mediated 
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Figure 6: Mass contours for xi (blue), X 2 (gray) and X 3 (red) in the (rriD, ms) plane for 
different choices of y and tand. In the green region all three Xi have a sizeable doublet 
component, see the text for details. 


by t-channel 'ip'^/Xn exchange. This behaviour is illustrated in the scatter plot of fig. 
where we can see that the avww is maximised for nearly degenerate spectra. 

Translating the bounds obtained in the previous subsection into bounds on 8 m = 
m^+ — we have: For y < 0.01, m^ < 280GeV and jf/np < 0.5 the exclusion is 
5m < 0.8GeV. For y = 0.2 the constrained region m^ < 220GeV, iVnp < 0.65 implies 
the exclusion of 5m < 20 GeV. For y = 1, since the exclusion of m^ < 275 GeV holds even 
for X being almost a pure singlet, the whole compressed regime is excluded. On the other 
hand, for positive tan0 and m^ < ITeV with 0.015 < \U\\\^ < 0.95, the lower bound on 
the singlet component of % translates into exclusion of 5m > 15 GeV. 


4 Light DM regime and LHC phenomenology 

In the previous sections we described the properties of the SDM and the constraints derived 
from DM direct and indirect detection experiments that, as we saw, can probe a very wide 
range of dark matter masses. In this section we would like to focus on the low DM mass 
region, as this is the relevant scenario for what concerns possible collider signatures. In the 
low mass region, besides the limits set by the DM searches, we implement other constraints 
on the model arising from collider experiments, in particular LEP and LHG. This makes 
manifest that the constraints on the model from DM searches and collider experiments are 
complementary, in particular in this regime. 
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In order to show the behaviour of the spectrum in the low-mass region of the parameter 
space, in fig.j^we plot the mass contours for Xi, 2,3 in the (tUd, rUg) plane. We do not display 
the region with < 100 GeV, because it is excluded by searches for charginos at LEP 
j38] . We also highlight in green the regions where all three states have a sizeable doublet 
component and thus can be relevant for the LHC production modes: >0.5, 

i = l,3. 

For small values of y (left panel of fig. |^, the second heavier neutral fermion is always 
mainly doublet, with a mass close to ttId, such as the charged states. The DM candidate xi 
is mainly singlet in the region with ttId >> mg so that ~ mg, while it is mainly doublet 
with ~ rriu for mg >> mo. The xs mass and mixing have a complementary behaviour. 
In the region rrij^ ~ rUg the three states have comparable mass, with X 2 mainly doublet 
and Xi and Xs largely mixed. These features do not have a sensitive dependence on tand. 
Such simple pattern is highly distorted for large values y, as shown in the right panel of 
fig-ii as a consequence of the large off-diagonal entries in the mass matrix. In particular, 
large mixing as well as large mass splittings among the three neutral states arise (and with 
respect to the charged ones). This might lead to interesting signatures at colliders, as we 
will discuss below. 

4.1 Dark Matter and EW constraints 

In this subsection we focus on possible modifications induced by the new physics sector to 
the measured properties of the SM particles. In particular, we quantify the contributions 
to EW precision observables (EWPO) and the invisible decay widths of the Z and h. 

• For what concerns EWPO, we use the expressions in isniiii] to compute the new 
physics contributions to the oblique parameters, AS and AT. As already observed 
in [12], the strongest constraints to the model arise from T, while the contribution 
to S is always small. Since AT scales as ~ {yf — y^Y, EWPO constraints can be 
relevant in scenarios with large y and large tan^. Given the small contribution to S', 
we simply impose the bound |Ar| <0.2 [iH] . 

• In the SDM, the Z boson width can receive relevant corrections if the decay Z XiXi 
in kinematically accessible. The contribution to the invisible width of the Z boson is 

= (15) 

where the Z’xiXi coupling is given in eq. pg. The Z width measurement performed 
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Figure 7: Constraints on SDM in the low-mass range for y = 0.2 and tan0 = 2, —20. 
Gray: region excluded by direct detection experiments. Purple: region excluded by indirect 
detection experiments. Blue: exclusion from BR(/j —)• invisible). 


by LEP impose the following constraint on possible new physics contribution to 
T{Z ^ invisible): AP^'' < 3 MeV (95% CL) jUj. 


• Analogously, the Higgs boson width into invisible particles can also receive important 
corrections in the SDM if < mh/‘i- The Higgs boson partial width to xi reads 


r(h ^ xixi) 


mh 

dyr 



( 16 ) 


where the coupling is as in eq. (10). New invisible decay modes of the Higgs would 
reduce all visible branching fractions, hence the invisible branching ratio can be 
constrained by fits to the observed decay rates. In the following we employ the 
following limit: BR(/j —)• invisible) < 26% (95% CL) ji2] . 


In the figures and we show the above constraints together with those from DM 
experiments discussed in the previous sections for the low-mass region of the parameter 
space. Regions excluded by direct detection are shown in gray, while the indirect detection 
bounds are collectively marked in purple. As before, the red line denotes = 0.12. 

From the figures, one can observe the resonant enhancement of the DM annihilation cross 
section when ~ mzl2 and ~ vnh/2^ respectively due to Z and h s-channel 
exchange. For these configurations, an efficient DM annihilation can be obtained even for 
a mostly singlet-like xi, i.e. mo ^ mg. 
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Figure 8: Same as in figure for y = 1 and tan0 = —20. Additional excluded regions 
that appear are the following. Green: exclusion from AF™'". Yellow: limit from EWPO. 


In figure we show the case y = 0.2 that resembles the Higgsino-Bino system of the 
MSSM, as discussed in the previous sections. Given the moderate value of the couplings, 
the EWPO do not pose any bound, whereas the Higgs invisible width constraint (blue- 
shaded area) can exclude regions with small in the case of low tan 6*. Larger and/or 
negative values of tan 6* can reduce or even evade this bound, as shown in the right panel 
of the figure. This is due to the dependence of the hxiXi coupling on tand that we can 
obtain from eq. (10). Furthermore, we see that there is a partial cancellation between 
the two terms of Ch^x case yi and y 2 have opposite signs. This is also responsible for 
the reduction of the direct detection bound in the right panel of the figure, as well as the 
decreased efficiency of DM annihilation on the h resonance [13] . 

The case of large y is strongly constrained by direct searches for dark matter. In 
particular, for moderate and/or positive values of tan 6, the LUX bound can be evaded only 
for ttId in the multi-TeV range. For this reason, we only show the case y = 1, tan 6* = —20, 
see figure As we can see, limits on invisible Z and h decays set relevant constraints for 
low values of rUg and nio, depicted respectively as green and blue regions. The bound from 
AT (shown as a yellow line) is also relevant although somewhat weaker. As mentioned 
above, DM experiments are very constraining, especially because of the increase of the 
couplings Chxx ^-iicl czxx large y. However, we stress that, as discussed in the previous 
section, the bounds obtained from DM experiments assume that 100% of the observed dark 
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matter is accounted for by our candidate xi- Thus they can be relaxed in the regions where 
a thermal Xi would give <C 0.12, assuming that it is only a subdominant component 
of the total DM abundance. In particular, in the special case ~ mh/2, where dark 
matter annihilation process is resonantly enhanced, the bonnds from direct and indirect 
searches can be evaded to large extent. For this reason, in the following we also discuss a 
distinctive collider signature that can be realised in such peculiar region of the parameter 
space (together with the more standard LHC phenomenology of the moderate y case). 

As a hnal remark, we can now complete the discussion of snbsection |3.2| setting com¬ 
prehensive bonnds also in the low mass region, considering at the same time all the ex¬ 
perimental constraints: Combining both analyses, we see that for y = 0.2, one can ex¬ 
clude X for 80GeV< < 220 GeV for Uh < 0.65. For the large Yukawa benchmark 
(r/ = 1), the bounds from invisible decays allow to keep the bound of the previous section 
{m^ < 275 GeV for Uii < 0.8) even for x of a few GeV. 


4.2 LHC signatures and constraints 

In this subsection we discuss the LHC phenomenology of the model. A detailed analysis 
and a recasting of existing LHC searches to this scenario are beyond the scope of this 
work. Here we restrict to a discussion of the generic features that are characteristic of 
different parameter space regions. An early study of the LHC signatures and prospects 
of the model can be found in [12] . As discussed in section the SDM can be considered 
as a generalisation of the Higgsino-Bino system of the MSSM. This analogy is very nsefnl 
in exploring the collider phenomenology of the model, since the snpersymmetric case has 
been extensively discussed in the literature with the 

important difference that here the mass splittings are controlled by y instead of g' and thus 
can be substantially larger than in the MSSM. 

The main production mechanisms at LHC are the electroweak Drell-Yan processes, 
i.e. charged fermion pair production throngh s-channel Z-boson or photon exchange, as¬ 
sociated production of charged and neutral fermions through s-channel W and neutral 
fermions pair production through Z\ 

pp , pp Xi'fp^ , PP XiXj , bi = 1,2,3 . (17) 

The production involving neutral states is controlled by their couplings to the W and Z 
bosons that, as we have seen, are weighted by their doublet components, i.e. by entries 
of the mixing matrix U, see eq. Q. In principle each of the three neutral fermions can 
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Figure 9: Heavy neutral states decays for y = 0.2 (left) and y = 1 (left) tand = —20. 
First row: contours of BR(x 2 ^ Xi^)- Second row: regions showing the dominant decays 
of X3- 


have a large doublet component, independently of the mass hierarchy among them, as 
has been highlighted in hgure The figure encodes the information about the dominant 
production modes we expect at LHC through the parameter space. Within the green 
bands all three neutral states can be substantially produced (obviously if their masses do 
not exceed the few-hundreds GeV range), while on the right (left) of green regions the main 
production modes involve the two heavier (lighter) states X2;X3 (Xi;X 2 )- The final state 
topology depends on the decay modes of the produced particles. The charged fermions 
will predominantly decay to Xi and a (possibly off-shell) bF-boson. This decay occurs 
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trough the (even tiny) doublet component of xi and it will be prompt in the mass range 
we consider, as soon as U 12 or U 13 > 0(10“®), i.e. y > 0(10“^). The decay modes of the 
mostly-doublet neutral fermions can instead present a rich pattern. The heaviest can have 
the following decay modes: 

X3 , X3 ^ hx 2 , X3 ^ X 2 Z , X3 ^ Xih , X3 ^ Xi-^ (18) 

Among these, the decay into and W is especially interesting. When this decay is 
kinematically open, the large conpling to W, c.f. eq. ([^, makes it the dominant mode 
leading to characteristic collider signatures, as we discuss below. The second heavier neutral 
fermion mainly decays to the lightest neutral fermion: 

X2 ^ Xi^ , X2 ^ Xi^ (19) 

and possibly to the charged fermion in very peculiar corners of the parameter space, where 
m ^2 > ‘^w + 'nT'D- Analytical formulae for the tree-level decay widths are reported in 
Appendix Obviously if the above modes are kinematically forbidden, the decays will 
occur through off-shell W, Z or h, depending on the conpling strength. 

In Figure]^ we show the dependence of the branching ratios on and rUg for r/ = 0.2 
(left) and y = 1 (right) and tand = —20. In the first row we plot contours of BR(x 2 —> 
XiZ)i in the second row we show regions corresponding to the dominant decay modes of 
X 3 . Gray-shaded areas highlight regions where only three-body decays are kinematically 
accessible, such that a reduction of the LHC sensitivity is expected. The blue-shaded re¬ 
gions are excluded by the collider constraints discussed in the previous subsection, namely 
EWPO, Z and Higgs boson invisible decay widths. 

The most interesting and easiest to probe final states at colliders are the ones involving 
the production of gauge bosons, that we will discuss in more detail in the following. Indeed, 
they can lead to distinctive and clean signatures at LHC, such as multi-leptons plus missing 
energy. We will focns here on the leptonic decays of the gange bosons and classify the 
possible signatures in terms of number, sign and flavours of leptons in the final statej^ 

®We consider only the cases corresponding to dominant combinations of productions and branching 
fractions. For instance, we neglect modes involving Xi X 2 -^, as Xi will typically dominate 

in the regions where both decays are kinematically open. Furthermore, for doublet-like X 2 and X 3 , the 
combination of mixing angles in czx^xn^ cf. eq. (|^, makes the production of X 3 X 2 always larger than XaXs 
and X 2 X 2 ) hence we only consider pp —X 3 X 2 among the neutral fermion production modes. 
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2 leptons. Events with two (possibly different-flavour) opposite-sign leptons plus missing 
energy can follow the charged fermion pair production 

pp —)• W^Xi^~Xi 1 (20) 

exactly like in the case of chargino production in the MSSM. Due to the large SM back¬ 
ground, this is one of the few electroweak production modes for which LHC experiments 
have not set stronger constraints than LEP yet. 

3 leptons, with a pair reconstructing the Z-boson mass. This final state leads to 
the most stringent LHC constraints at the moment. The production modes contributing 
to this signature are 

PP X2,3'tp^ ZxiW^Xi ■ (21) 

Hence schematically, the W plus Z plus missing energy production rate will be given by 

X BR(x3 Zxi) + criX2'4’^) X BR(x 2 ^ Zxi) . (22) 

Searches based on this kind of events have been published in |53l [5l] . The ATLAS search 
[S3] has been recently recasted in the context of the Higgsino-Bino system of the MSSM 
[33|- We can use this reinterpretation to estimate the LHC sensitivity to our parameter 
space. Erom the analysis of (33], we see that the ATLAS search [S3] is sensitive to our 
model roughly under the following conditions: when < 270 GeV, ^75 

GeV and the average of BR(x 2 ^ -Z^Xi) and BR(x 3 —>■ Zxi) is at least 60%. Let’s consider 
for instance the case y = 0.2, tan0 = — 20 . From the mass contours and the branching 
ratios shown in figures and we see that the region relevant for the Z and h resonances 
in figure with m-o < 270 GeV is likely excluded by the ATLAS search [53|. For the 
other cases shown in figure and the bound coming from BR(h —>■ XiXi) seems to be 
at present more stringent than the 3-lepton search. 

3 leptons, without a pair reconstructing the Z-boson mass. This final state can 
be obtained from decays of three W bosons. The production mechanism is X3'^^ with X 3 
decaying to 


pp ^ W^W^XiW^Xi ■ (23) 
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This case can be realised if > mw + n^D and xs has a sizeable doublet component. 
These conditions are difficult to obtain for low y, as we can see comparing the plots of 
figures 1^ and in the region where the decay of X 3 into W is preferred, the doublet 
component of X 3 is small, reducing the production cross section. Instead, this signature 
can be realised in the case of large y. For instance, consider the region of the plane in 
figure]^ where ~ vah/^ in the case oiy = 1 and tan0 = —20, e.g. with rUg — ~ 150 

GeV. There Xs is mostly doublet and decays preferably to realising the scenario 

we just described. 

4 leptons with one pair reconstructing the Z-boson mass. This final state can be 
obtained via the process: 

PP X 3 X 2 W^'ip'^Zxi W^W"^xiZxi (24) 

This case can also be realised in the same region of the parameter space of the y = 1 and 
tan 6 * = —20 case discussed above, just close to the line of Higgs enhanced annihilation in 
figureIndeed there X 2 is mostly doublet and decays preferably to Z plus Xi, see figure]^ 
The final state signature is two W bosons (of opposite sign) plus one Z and missing energy. 
Two out of the four leptons will have same flavour and opposite charge and reconstruct the 
Z boson mass, while the other two will have opposite charge but of uncorrelated flavour. 

4 leptons with two pairs reconstructing the Z-boson mass. This case follows from 
the same production mode as the previous one when xs also undergoes the more ordinary 
decay to Z: 

PP X 3 X 2 ZxiZxi ■ (25) 

Considering the behaviour of the branching ratios shown in figure this mode is particu¬ 
larly relevant in the regions with small rUg and rUn in the few hundred GeV range. This 
makes it, together with the 3-leptons mode from WZ discussed above, very promising to 
test the unconstrained regions in figure and at the current LHC run. 

5 Summary 

In the present work we have analysed the theory and the experimental constraints of a 
dark matter candidate that is a mixture of an electroweak doublet and singlet fermion. The 
model enjoys certain properties: i) It represents a minimal viable description of electroweak 
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dark matter that is stable due to a symmetry of the theory, ii) The currently running 
direct detection experiments are sensitive enough to probe DM-nucleon interactions via 
Higgs exchange. Since this model has tree level DM-Higgs interactions, it can be used 
as a simplified model to parametrise the bounds on this coupling, iii) It improves gauge 
coupling unification of the SM. iv) It captures the relevant dynamics of DM candidates of 
UV theories such as the Bino-Higgsino system in MSSM and the Singlino-Higgsino system 
in NMSSM-type scenarios. 

We have used a combination of results from direct search experiments (LUX [2^, PICO 
|27jb indirect search experiments (FERMI [26], ICECUBE [28]) and collider searches 
(Higgs and Z invisible decays as well as EW precision measurements) in order to set 
constraints on the viable parameter space of this model. 

We focused on three different regimes for the DM-Higgs Yukawa coupling y. Small 
(y = 0.01), MSSM-type {y = 0.2) and large {y = 1). Lowering y below 0.01 does does not 
change substantially the phenomenology of the model while values relatively larger than 1 
are not compatible with perturbativity up to the unification scale. 

We have seen that, while colliders can probe DM candidates up to a few hundred GeV, 
direct and indirect searches allow us to extend the limits up to masses of several TeV, 
depending on the composition of DM. The bounds on the parameter space from direct and 
indirect DM searches are summarised in figures [T] -while the combined bounds from all 
experiments in the low mass regime are shown in figures and For small y (< 0.1) the 
FERMI limits on %% —?• WW exclude % that is mostly doublet up to ~ 280 GeV, while 
for y = 0.2, SI, SD and collider limits are also added and allow one to exclude the DM 
candidate y with a singlet component < 0.65 for masses < 230 GeV. Eor y = 1, SI 
direct searches set strong bounds on the model, apart from regions of the parameter space 
where the DM-Higgs coupling is vanishing. We have shown that these ‘blind-spots’ are 
probed by other searches, in particular SD and indirect searches. This allows to exclude y 
with mass < 275 GeV unless it is a highly pure singlet state {Uh > 0.8). If we restrict to 
positive values for the Yukawa couplings, y is excluded up to several TeV unless it is an 
almost pure state. 

We have also discussed the phenomenology of the model at LHG. The signatures are 
similar to the more studied Bino-Higgsino system with the important difference that the 
coupling to the Higgs and the mass splittings are controlled by y instead of the hypercharge 
coupling. The most interesting final states involve 2, 3 or 4 leptons plus missing energy. It 
would be interesting to perform a detailed study of such signatures at LHC at ^/s = 13 TeV 
and in particular to ascertain the reach for the region of parameter space that is still allowed 
in the moderate-low mass case (see the white region of figure . 
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Regarding future directions, in this work we have not included constraints coming from 
limits on fluxes of charged cosmic ray particles. It would be interesting to see how these 
limits can be used in combination with the present ones in order to further constrain the 
model. Also, the computation of DM-nucleon scattering amplitude has been done at tree 
level. A more detailed treatment would require the inclusion of loop corrections that can 
be relevant [SH] in regions of the parameter space where the tree level coupling is small, 
such as the pure doublet limit. 
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A Decay rates 


In this appendix we report the analytic formula [HS] for the decay of the charged fermion 
and the heavy neutral fermions used in the main body of the paper. 

Given a generic Yukawa-type interaction of the form 'ipiciPL + CRPR)x(t>, the partial 
width for the decay % —?• is given by 


r(x ^ 


327rm^ 


y/2 


{{\cl? + \cR\^){m\ + - ml) + + cl 4 )) , 

(26) 


where A(a, b, c) = + b‘^ + — 2ab — 2ac — 2bc. 

Given a generic vector boson interaction gV^yx^iyCRPi + crPr)x the partial width for 
the decay y —>■ is given by 


Tix^vy) = 




I,ml,ml 


)l/2 


327rm^ 

X 


-6 {crcI + CLcy)m^m^+ 


(|clP + |ci^p)(m^ + ml-ml + {{ml - mlf - ml)/ml)] . (27) 


We are interested in the possible decay of the mass eigenstate fermions in the singlet- 
doublet model described in Section The coupling of the charged fermion with the three 
neutral fermions and the W boson is given by 


where here i refers to the i-th mass eigenstate of the neutral fermions, ordered in absolute 
value of the mass. The couplings among the neutral fermions and the Z-bosons are 


(cl)., = -^(UhU,^ - Ul,U,s) , (cr)« = -(cl ),,. 


(29) 


Finally, the coupling to the Higgs boson are[^ 


(cl)., = -^ {viUiu;, + y,u;,ui + y^uiu;, + v,uiu;^ , (ci,., = ( 4 )« ■ (30) 


^Compared to the expression in the main body of the paper, we have symmetrized and taken into 
account a 1/2 factor in the Lagrangian. 
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